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STEADY NUCLEAR COMBUSTION IN ROCKETS* 
By E . Sanger 


SUMMARY 

The astrophysical theory of stationary nuclear reactions in stars 
is applied to the conditions that would he met in the practical engi- 
neering cases that would differ from the former, particularly with 
respect to the much lower combustion pressures, dimensions of the react- 
ing volume, and burnup times. 

This application yields maximum rates of heat production per unit 
volume of reacting gas occurring at about 10® °K in the cases of reac- 
tions between the hydrogen isotopes, but yields higher rates for heavier 
atoms. For the former, with chamber pressures of the order of 100 atmos- 
pheres, the energy production for nuclear combustion reaches values of 

about 10 4 kilcalories meter‘d second"- 1- , which approaches the magnitude 
for the familiar chemical fuels. The values are substantially lower for 
heavier atoms, and increase with the square of the combustion pressure. 

The half-life of the burnup in the fastest reactions may drop to values 
as low as those for chemical fuels so that, despite the high temperature, 
the radiated energy can remain smaller than the energy produced, partic- 
ularly if an inefficiently radiating (i.e., easily completely ionized 
reacting material like hydrogen), is used. 

On the other hand, the fraction of completely ionized particles in 
the gases undergoing nuclear combustion must not exceed a certain upper 

limit because the densities (~10" 1( ^ g cm"^) lie in the range of high 
vacua and only for the previously mentioned fraction of nonionized par- 
ticles can mean free paths be retained small enough so that the chamber 
diameters of several dozen meters will suffice. 

Under these conditions it appears that continuously maintained stable 
nuclear reactions at practical pressures and dimensions are fundamentally 
possible and their application can be visualized as energy sources for 
power plants and propulsion units. 


"Stationare Kernverbrennung in Raketen." Astronautics Acta, I, 
Fasc. 2, 1955, pp. 61-88. 
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Three basic methods of utilizing the enthalpy of nuclear-reaction 
gases, of the order of 10^® kilocalories per kilogram, in jet-propulsion 
units, suggest themselves: (l) the direct expansion of these gases in 

a pure atomic rocket; (2) the admixture of the reacting gases with the 
surrounding air in turbojets, ram jets, or with other inert gases in 
thermal atomic rockets; and (3) the conversion of energy in photon gas 
in the photon rocket. 

For the last two methods, the combustion pressures available in 
practice seem insufficient to attain the required chamber loading. 

On the other hand, the technical conditions for the operation of 
stable nuclear reactions in central power stations , marine propulsion 
plants, or in pure atomic rockets of very large dimensions seem much 
more attainable. 


I. SOCIOLOGICAL AND TECHNICAL ASPECTS OF THE PROBLEM 

Technical realization of space travel demands, on the part of 
humanity, an effort that may be represented by an eleven place figure 
of working hours. In the course of a lifetime, an individual can at 
best contribute a five place figure; therefore, several million people 
will have to occupy themselves principally with space travel. This is 
already the case if air travel is regarded as an initial stage of space 
travel . 

This effort towards space travel can be made in two ways: 

(l) either by a vast, concentrated program; or (2) , by many small con- 
tributions. These programs are similar except that the second effort 
presents different psychological problems. 

Much strenuous effort would be required to achieve the technical 
development of manned space machines on the basis of presently existing 
chemical rocket motors. This method, while still in the initial stages 
has been followed under the pressure of military necessity. However, 
majority opinion hesitates at so concentrated an undertaking, since the 
investment is so high for so limited a goal and since the millions of 
people needed to work on space travel must include not only the pro- 
fessional people working on air travel but also those who could better 
serve in other occupational categories. 

The less spectacular path of many smaller contributions to technical 
research, especially in the field of nuclear- powered rocket motors, prom- 
ises to be less expensive, less risky, and more effective than the vast, 
concentrated program. In this research, the number of people required 
would be limited to those now working in air travel and in similar 
fields . 
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This second method of achieving space travel has already been fol- 
lowed, since a major aspect of the problem (i.e., the method of obtaining 
propulsive energy from the conversion of mass to energy) is the present 
concern of the people working with power development, weapons, land and 
marine transportation. 

Both methods are in a healthy competition with each other, but the 
former, which makes use of chemical reactions, attains a conversion of 

but 10" ^ of the mass into energy, whereas the second method strives for 
the 100-percent conversion occuring in nature but presently achieves a 

conversion of about 10”^ in nuclear combustion. When technical mastery 
is achieved, this method of producing energy promises greater repercus- 
sions on the history of humanity than were produced by the deed of the 
mythical Prometheus. 

According to the views of the astrophysicists, steady nuclear com- 
bustions are the predominant source of energy release in the universe 
and thus the source for all life on earth. The state of nuclear physical 
research permits us to assume that the new Prometheus already lives among 
us who will also steal fire from heaven and bestow it upon humanity for 
immediate use as steady nuclear combustion. 

The difficulties of technical nuclear combustion arose primarily 
from the fact that combustion pressures and burnup times must be much 
smaller than the values attained in the cosmic burnup at the centers of 
the sun and stars. It is with these technical difficulties, in partic- 
ular in the application of the steady nuclear combustion in rocket 
motors, with which we will be concerned. 


II. THE NUCLEAR COMBUSTION PLASMA 

A gas at rest contains, according to the Maxwell-Boltzmann statis- 
tical distribution, a most probable thermal particle velocity 
C 2 = 2kT/AM, determined by its temperature T with a distribution of ther- 
mal velocities of the particles extending from zero and infinity. 

The collisions between the gas particles, under the assumption of 
normal temperature of the gases, will be predominantly of an elastic 
nature (i.e., there is neither loss nor gain of translational kinetic 
energy) . 

If, within the Maxwell distribution of the given gas we pass to 
higher velocities, we soon encounter inelastic impacts in which energy 
is lost by excitation of internally quantized degrees of freedom, namely, 
rotation, vibration, or excitation of the gas molecules. Energy is 
thereby removed so that the translational energy and the translational 
temperature of the gas particles is reduced. 
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This energy transfer to the internal degrees of freedom can, in 
turn, give rise to the emission of photons which, if the gas volumes are 
assumed to be infinitely large compared with the photon mean free paths, 
form a photon gas in equilibrium with the corpuscles of the rest mass that 
constitutes the gas itself. 

As we continue to pass to the more energetic but less frequent 
collisions, we find the dissociation of gas molecules into radicals and 
atoms. These in subsequent inelastic collisions, with suitable recom- 
bination partners, may loose their internal energy so that the transla- 
tional energy and, therefore, the temperature of the gas, increased. 

With still higher collision velocities complete ionization of the 
atoms is possible, which then opens the way for collisions between the 
nuclei themselves. 

These nuclear collisions will at first again be predominantly 
elastic. Only in very rare cases the still relatively small collision 
velocities will lead, as a consequence of the wave mechanical tunnel 
effect, to nuclear reactions. These collisions will yield on the aver- 
age up to 10^ times higher energy release than the chemical recombina- 
tion collisions. Only for those very rapid and, therefore, extraordi- 
narily rare collisions whose energy is sufficient to overcome the clas- 
sical Coulomb potential of the nucleus, will there be an appreciable 
probability for nuclear collisions that lead to nuclear reactions. 

These reactions will occur continuously in the same way as the chemical 
reactions of combustion gases, but will occur so rarely as to be negli- 
gible from the energy production standpoint. 

With true adiabatic isolation of the gas under consideration, these 
very rare nuclear reactions must naturally also lead to a gradual in- 
crease in the temperature in accordance with the Maxwell distribution, 
which then leads to a more and more rapid increase in the frequency of 
all the processes described. Since the process accelerates at high 
collision velocities, an avalanche- like increase of the number of ther- 
mal nuclear reactions in the gas results, which tends to reach explo- 
sively an end temperature of the order of magnitude of lO^® degrees, if 
it consists of light elements. 

This natural and self-supporting process does not occur in our more 
restricted environment because no gas masses of infinite volume exist 
and no boundary of a finite volume is completely impervious to heat. 

If man causes this process artificially, then he will have to form 
boundary conditions to have as low a heat loss as possible and, at the 
same time, accelerate the process. This can be accomplished, for 
example, by reducing the long start-up period through artificial heating. 
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This method has already succeeded in bringing about explosion- like 
thermal nuclear combustions of certain kinds of atoms to a technically 
useful scale. 

The original form in which the energy appears in chemical reactions, 
and the mechanism according to which it is distributed over the combus- 
tion gas particles are as yet only very little known, and has only 
recently begun to be investigated. 

In the case of nuclear chemical reactions, we are far better in- 
formed on the primary heating up process itself owing to the direct pos- 
sibility of the observation of the paths and velocities of the nuclear 
reaction products as a result of their ionization action in the sur- 
rounding gases, a circumstance that is lacking in the case of the chemi- 
cal reaction products because of their much smaller energies. 

The nuclear chemical heat of reaction can readily be computed in 
advance, with sufficient accuracy, from the mass defect of the particles 
taking part in the reaction. This is impossible in the case of chemical 
reactions since the mass defects are too small to be weighed accurately. 

Furthermore, the multiplicity of the degrees of freedom for receiv- 
ing this released energy, in the case of chemical reactions (translation 
of the reaction products and possibly of photons, rotation, oscillation, 
different energy level excitation, and ionization), is considerably 
greater than for nuclear reactions where only the translation of the 
reaction products and of the photons are principally being considered. 

The nuclear reaction proceeds in such a manner that the two reac- 
tion partners combine to form a compound nucleus of very short life, 
after which, depending on the degree of excitation and their individual 
properties, they divide approximately in half (fission), into very many 
smaller parts (spalling), or emit individual small particles (a-particles, 
protons, neutrons, electrons, ppsitrons, and photons). The last process 
has been principally investigated in nuclear combustions at the present 
time. 

The excess energy (heat of reaction) appears here as kinetic energy 
of translation, the newly formed particles disperse with great velocity, 
as we also have recently been more and more forced to assume for the 
case of chemical combustion reactions. 

Since the nuclear chemical burnup process for readily attained 
pressures often lasts considerably longer than that of the chemical com- 
bustions of equal pressures, the newly formed particles experience very 
many elastic collisions with the remaining particles of the reaction 
plasma (i.e., with atomic nuclei and electrons), and with these particles 
the Maxwell-Bolt zmann velocity distribution is immediately assumed. 
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For the temperatures and pressures that are of interest here, the 
mean free path of the particles is of the order of magnitude of 10 7 
centimeters within a fully ionized plasma and 10" 1 centimeters in a non- 
fully ionized plasma of the same temperature. In a surrounding cooler 
gas jacket in which nuclear reactions occur more rarely, the mean free 
path is many tenth powers less. 

The particles, therefore, will instantly diffuse from a fully 
ionized plasma of steady pressure, even if the diameter of the combustion 
zone was of astronomical dimensions. 

Only in incompletely ionized plasmas obtained (e.g., by the inclu- 
sion of heavy ions) , can the dimensions of the reactions remain feasible 
even at the highest possible pressures. 


III. THERMAL NUCLEAR REACTION AND ENERGY RATE 

The problems of the production of energy through thermal nuclear 
reactions have been clarified particularly through the work of the astro- 
physicists and we here follow essentially the known work of G. Gamow 
(ref. l). 

Nuclear reactions between fully ionized atoms (i.e., bare positively 
charged atomic nuclei), become very probable when the relative collision 
velocities of the particles are so high that the centers of the nuclei, 
during the collision, approach against the repulsion of the Coulomb field 
to a distance equal to the sum of the de Broglie wavelengths of the 
particles. 

From considerations of wave mechanics the potential barrier in some 
cases can be penetrated by slower collisions, although with a probability 
which decreases exponentially with the collision velocity. The proba- 
bility of the nuclear reaction by a particle which has penetrated into a 
nucleus is, in addition, determined by the resonance relations of the 
compound nucleus formed, and by the escape probability of the particle 
from this intermediate nucleus. 

All these influences on the reaction probability are evident in the 
behavior of the so-called reaction cross-section a, the variation of 
which is plotted schematically in figure 1, against the collision veloc- 
ity, for reactions between charged particles. 

The cross section for thermal nuclear reactions between charged 
particles increases exponentially with the collision velocity w until 
the kinetic energy is of the same order of magnitude as the potential 
energy at the summit of the potential barrier; the increase then becomes 
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slower, the cross section reaches a maximum, and then,, at very high 
collision velocities, it decreases with further increasing v because 
the duration of time that the colliding particles remain in the hit 
nucleus becomes increasingly shorter and less sufficient for a nuclear 
reaction. 


Gamow gives the following quantitative expression for this relation: 


cr JPjt exp(8^eV2RAMZ^Z^ |h - rc 2 e 2 Z 1 Z 2 |hw) 4jt^R 2 AM|h 2 [cm 2 ] (l) 

for which the following symbols are used: 


X = \/2n 
X « h/AMw 

A a Aikz/(Ai + A 2 ) 

A 1> A 2 

Ma 1.6604x10" 24 g 

e = 4 . 805 xlO -10 cm 3 / 2 g 1 / 2 sec" 1 

R a 1.7 + 1.22(A X + A 2 ) 1 / 3 X10" 13 cm 

Z 1> Z 2 

h a 6.624X10" 2 "'' erg-sec 


scattering wavelength, cm 

de Broglie wavelength of material 
wave of particles, cm 

reduced atomic weights for energy 
absorption in plastic collision 

[-] 

atomic weights of collision 
partners [ - ] 

nuc leone mass 

elementary charge 

radius of compound nucleus 

atomic numbers of colliding 
particles [-] 

Planck's constant 


T half width of nuclear resonance 

level [erg] 

The three multiplied factors in equation (l) have the following 
physical meaning: 

denotes the effective scattering cross section; the exponen- 
tial represents the transfer coefficient, which also contains the tunnel 

effect. For the collision velocities greater than w’ = eit/2 -y/Z-^Zg/RAM 

the exponential is greater than one, so that the reaction cross section 
a can become fundamentally greater than the scattering cross section. 



8 


NACA TM 1405 


The maximum value of the transfer coefficient occurs for w = 00 and is 

of the order of magnitude 10 4 . The velocity w* is equal to that which 
is required for the nuclear centers to approach to a distance 8/rt^R *»R. 

Finally, the fraction appearing as the third factor is the ratio 
of the reaction probability 2«r/h to the eigen frequency h/2 rtR 2 AM 
of the nucleus. This ratio for (p,a) reactions lies in the neighborhood 
of one, so that each penetration of the potential wall leads to the 
reaction. For (p,y) reactions, the ratio is generally 10 4 to 10^ times 
smaller, and for (p,P) reactions, up to 10 24 times smaller. Only by 
experiment in certain cases can f be more accurately determined. 


Equation (l) can be reduced to a somewhat more manageable form: 


a « R jt 


o p Zp ”1 p 

^ AM~2 6XP h — ~ (2-VRAM - rte-Zz^g/w) [cm* 1 ] (la) 


The variation of the reaction cross section o with the collision 
velocity w shows a maximum at the collision velocity 


( w )cr v “ 2 V^ n 2 e 2 Z 1 Z 2 /h [cm s" 1 ] 


which, thus, has no connection with the velocity w* mentioned pre- 
viously. The maximum cr corresponding to this velocity is 


°max “ 


Sn^AMe^Z^Zg 


/8rte-y/2RAMZ 1 Z2 \ 


The optimum collision velocity with respect to o according to 
equation (lb) can be obtained from the following considerations. 


The absorbed kinetic energy corresponding to (w ) n is given by 

u max 


AM(w 2 ) 


4* 4 AMe 4 Z 2 Z 2 


and the corresponding de Broglie wavelength: 


y. = — = h/2« AM(w) 


4-^2 n 5 AMe 2 Z 1 Z, 


Therefore, 


Z n Z 0 e< 


V2 /** 
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The absorbed kinetic energy is equal to the Coulomb potential energy 
at the mutual approach of the particles to a distance of -y/2/jt times 
the de Broglie scattered wavelength. 

If the cross section variation in figure 1 is known, the number of 
nuclear reactions taking place in a given plasma may be computed by 
multiplying the probability curve of the reactions (a) by the probability 
curve of the thermal velocities ( dN / dw) . The product indicated as a 
dotted curve in figure 1 has a sharp maximum at 


w opt 






AMh 


( 2 ) 


Gamow integrated the product versus w by replacing the actual 
curve with a Gaussian probability curve of the same height and width and 
thus obtained^ for the energy production through thermal nuclear reac- 
tions per unit mass and unit time of the reaction plasma 


where 


AE — Epc 1 c 2 h 1 


(k 2 T - 1 / 3 ) 2 
exp^k 2 T -1 / 5 j 


, erg g sec 


-1 


( 3 ) 


E 

P 


heat of a single reaction; ergs 
density of plasma; g cnT^ 


c l' c 2 


relative weight concentrations of the reac- 
tion particles 


k — 1.3805 *10" erg-deg - ^ Boltzmann constant 


2 ThR 

k x = -575 “3 2 exp 81 

3 d / c AAt Ac>e^Zi Z 


U A 2 e *1^2 


f2rt 2 RAMe 2 Z 1 Z 2 \ 1 / 2 

V ? ) ’ 


4„2„2 


ko s 54n' 


j Me ZpZ 2 




deg' 


1/3 


cm^ g“2 sec -1 (3a) 


(3b) 


kh 


k-^ and k^ are the constants determined by the nature of the reacting 
particles and may be computed once and for all. 
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Differentiating equation (3), there is also obtained a maximum of 
the energy rate with respect to the temperature at 


, 54jt 4 AMe 4 Z?z| 
T opt - (V 2 ) - - 


(3c) 


8kh 


and of an amount 


AEmax 15 ( 2 / e ) 2 E P c i c 2 k l 


8 _ „ R 2 rh 

-57F EC 1 C 2P 


it 


M 3 AA 1 A 2 e 2 Z 1 Z 2 


(' 

QXpl- 


QjtG-y^ 2RAM 


2 ) 


(3d) 


The position of this optimum, which is independent of the density 

of the plasma, is approximately obtained from the consideration that in 

figure 1 the position of (w)^ , according to equation (lb), is inde- 

max 

pendent of the temperature, whereas the most probable velocity c in- 
creases with increasing temperature T and can be made to coincide with 

(w) . If, therefore, the plasma temperature is so chosen that the 

u max 

velocity corresponding to the maximum of the reaction rate coincides 

with (Azc. /A.)' 1 / 2 times the most frequent velocity c = (SkTz^/A^M) 1 / 2 
i 1 1 i 

of the Maxwell distribution equation (lb) , then this is approximately 
the temperature of the greatest reaction probability, thus 


4jt 4 AMe 4 Z 2 z| 

^ S 2 


(4) 


Equation (3c) gives a more accurate value, about 69 percent higher, 
arising from the fact that the area under the adN/dw curve does not 
accurately become the greatest when the maximum of this curve is the 
greatest. 

If the maximum of the Maxwell curve dN/dw in figure 1 coincides 
with the maximum of the a curve, then the adN/dw curve must also 
have its maximum at the same point. 


Finally, the very important optimum temperature of the thermo- 
nuclear combustion at constant density can, according to equation (3c), 
also be understood clearly from the following considerations . 


; tt^ 
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In the inelastic collision between two particles of different 
weights the absorbed energy, according to the classical laws of impact, 
is obtained from the momentum and energy laws as 

AiA2(Mw2) _ 

A x + Ag 2 2 V 1 

In order to make two charged particles approach each other in oppo- 
sition to the Coulomb repulsion to within a distance of Q/n-yJbl times 
the de Broglie scattering wavelength \ * \/2rt = h/2jt AMw, there is an 
accompanying absorption of kinetic energy of an amount 

Z^Z^^^ -^/siitZ-^Z^e ^ 2^:AMw AMw^ 

8/*V54xX Sh ~ 

and, therefore, the collison velocity 

-v/54« 2 e2z l Z 2 
w = E — 

If this required collision velocity w becomes equal to the most 
frequent velocity c Z = 2kT/AM in the Maxwell-Boltzmann distribution, 
then the gas temperature is again T 0 p-t, according to equation (3c) . 

In a manner similar to the method used to find the position of the 
maximum of the energy, production curve equation (3) we may try to 
estimate the value of this maximum. 

According to Maxwell-Boltzmann statistics, the number of collisions 
between different particles in a mixture of two kinds of gas per unit of 
mass and time is 


N « 


4c2C2pV^ 
M 5 / 2 A !A 2 VS5 


OS 


( 6 ) 


hence, at the temperature of the peak of the energy production curve, 
that is, with the substitution of equation (3c) 


-J 54ic 3 / 2 ciC2pe 2 Z]_Z2 

N - 5 0S 

AiAgM^h 


(7) 


If the greater number of the collisions do not lead to scattering 
but lead to a reaction, there is substituted in place of the scattering 
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cross section as the reaction cross section ~a_„„ as given by equa- 

LILcLa. 

tion (lc). To obtain the approximate value for the maximum of the energy 
rate curve. 


^max ' E= 1 c 2 P 


R 2 rh 

AA 1 A 2 M 3 e 2 Z 1 Z 2 


A 

expl- 


8ice- > /2RAMZ 1 Z 2 


- 2 ) 


( 8 ) 


which is thus about 2 percent less than the accurate value according to 
equation (3d) . 


With the previous considerations it is now simple to investigate 
the technically more important energy rate of the nuclear combustion at 
constant pressure p instead of, as up to the present time, at constant 
density p. 

If, using the gas equation for an expression of the density in 
terms of the 


P _ JL t 
p AM 


or 


pressure (where the Loschmidt number 
inated with the aid of the relation 
made in equation (3), the expression 
becomes 


pM/kT 

p “ 


L a 6.023 X 
LM « 1 g mol 
for the rate 


AE* = Epcj_c 2 


Mk n 


(kgT -1 / 5 )' 


( 9 ) 


10 mol was elim- 
a substitution is 
of energy production 


kZ c i/ A^ t exp(k 2 T“' L /^) 


erg g"- 1 - sec 


.-1 


(10) 


The energy production per unit mass and time at constant pressure 
p has a maximum with respect to the temperature T which is independent 
of the pressure at a now lower temperature: 

T 0 pt‘ * ( k 2 / 5)3 (l ° a) 

and of the amount 

^max' - ( 5 / e ) 5 Ep Cl c 2 1 — ( 10b ) 

kk 2? c i/^i 


I TT^ 



; tt?7 


NACA TM 1405 


13 




The length of the nuclear thermal burnup is characterized by the 
half value time t of the reaction: 


t = 


^max' 


In 2 


( 11 ) 


The energy production per unit volume at a given combustion pres- 
sure is of greatest practical interest in achieving nuclear combustion. 
This value is obtained from equation (3) through multiplication by p 
and replacing the density p by p and T by means of the gas equa- 
tion (9) to yield 


AE" 


M^kn 


Ep^c^c 2 


(kgT- 1 / 3 ) 


(kZcj/Ai) 


2 T 2 exp(k 2 T _1 / 3 ) 


erg cm sec x 


(12) 


This energy production per unit of volume and time for constant 
pressure of combustion p has, with respect to the temperature T, a 
maximum at the now still smaller temperature and is again independent of 
the pressure: 


T opt" - (V 8 ) 3 < 12a ) 

and of the amount 

o 9 M 2 k x 

^ax” “ (%) E P c l c 2 6 o 

k^kZc./A.^ 
i 

In table I, for fifteen chosen reactions between light, charged 
nuclei, characteristic numerical magnitudes of stoichiometric nuclear 
combustion computed from equations (l), (3a), (3b), (3c), (3d), (lOa), 
(ll), (12a) , and (l2b) are given. 

It should be emphasized now, however, how extraordinarily arbitrary 
this choice is relative to the thousands of possible two-particle impact 
reactions between the elements of the entire periodic system, how many 
times greater the possibilities are for three-particle reactions or for 
a still greater number of particles, and how the recently developed 
field of meson physics makes our present ideas appear as mere childish 
efforts . 

The first column of table I gives an index number to which the re- 
action is referred on later diagrams. 


(12b) 
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The second column defines the corresponding reaction in the usual 
notation. It is seen that several of the first elements of the periodic 
system are combined with protons, deuterons, or tritons. Practically, 
it is to be noted here, that in a nuclear combustion hardly a single one 
of these simple reactions takes place by itself, but there is generally 
an entire series of reactions, as for example, the familiar astrophysical 
carbon-nitrogen cycle of bright stars leading to the formation of helium 
from hydrogen (fig. 2(a)), the hydrogen chain of the darker stars simi- 
larly resulting in helium from hydrogen (fig. 2(b)), or the surmised 
hydrogen bomb process for forming helium from lithium hydride (fig. 2(c)). 

The third column shows the heat liberation of the simple reaction 
twice, first in the usual unit of the nuclear physicists in MeV per 
single reaction, and then in the engineering unit calories per gram (i.e., 
per unit mass of the propulsive material). In the latter figure it is 
immediately evident that the energy releases are higher by six to seven 
tenth powers than the chemical heat release values. 

The fourth column contains the measured half widths r of the 
nuclear resonance in eV. For the numbers put in parentheses no measured 
values are known and have been estimated partially on the basis of the 
considerations leading to equation (l). 

The fifth and sixth columns show the material constants of the Gamow 
equation for the energy rate according to equations (3a) and (3b) . 

The first results of technical significance start with columns 7 
and 8 where the maximal energy rate per unit mass AE^^/p at constant 

plasma density p and the corresponding optimum temperature are 

given. Since the energy rate is proportional to the density it is possi- 
ble to obtain the value for an arbitrary density from the table. The 
numbers shown represent the rate at the density p * 1 gram per cubic 
centimeter, which is close to the conditions as they may occur at the 
centers of the stars and in atom bombs. The half -value times of the re- 
actions, which are not given in this particular instance, go as low 

as 10" ^ second; the burnup under these conditions occurs in many re- 
actions extremely rapidly and explosively; the optimum temperatures are 
always essentially higher than those at the centers of normal stars. 

For continuous technical nuclear combustions, the succeeding columns 
are more important. Columns 9, 10, and 11 show the maximum energy rates 
at constant plasma pressure and again per unit mass, per unit pressure 
in atmospheres, so that the energy rates for any steady pressures can be 
directly obtained by multiplying the numerical values for the unit pres- 
sure by the corresponding pressure. The corresponding optimum tempera- 
tures are always lower here by more than a tenth power. 
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The half-value times for 1 atmosphere pressure are given; for other 
pressures they are obtained by dividing the given figures by the desired 
pressure. At the technically moderate pressure of 1 atmosphere, the 
half-value times of the reactions lie between a million times the age of 
our universe and fractions of a second, so that the choice of reactions 
for stationary technical nuclear combustions must primarily be made 
according to this column and not according to the heat of the reaction. 

Finally, in columns 12 and 13, the maximum energy rates per unit 
volume are given at constant plasma pressure, again for the unit pres- 
sure 1 atmosphere. At other pressures, therefore, the energy rate per 
cubic centimeter is obtained through multiplication by the square of the 
corresponding pressure. 


The pres sure -independent value of the optimum temperature is 
observed to be noticeably smaller compared with the values of column 9. 


The fastest reactions under these conditions are those of deuterium, 
tritium, and lithium at temperatures of the order of magnitude of 10® 
degrees for the lighter elements and 10® degrees for the heavier elements. 


The energy production per unit volume remains moderate in all cases 
because of the extremely small gas density. Even at a 100 atmosphere 
pressure, the energy rate barely reaches the value of 10^ kilocalories 
meter - 3 second* 1 for the economically practical nuclear materials and 
about 10 times more for the tritium combustion, as compared with the 
known combustion chamber loadings up to 10® kilocalories meter"^ second - ^ 
of chemical rockets. 


Here we encounter the first fundamental difficulty of continuous 
thermal nuclear combustion, the small burnup velocity and combustion 
chamber loading. 

So far we have restricted ourselves 'to the quantitative energy pro- 
ductions at the corresponding optimum temperatures and will now investi- 
gate the dependence on the temperature numerically. 

In figure 3 the energy rates per unit density of the nuclear reac- 
tions of table I are plotted against the plasma temperature T accord- 
ing to equation (3) . As stated previously, the energy rate per unit 
mass, because of the dependence on the number of collisions per second 
upon the mean free path, is proportional to the density. 

At constant density, the most favorable temperatures for energy pro- 
duction lie, as shown in table I, between about 10 9 and 10 12 degrees in 
all cases, increasing with the atomic number of the reacting particles, 
and independent of the plasma pressure. 
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The highest point of the production curves is essentially determined 
by k-L or by the maximum reaction cross section o max as given by 

equation (lc), and by E,p,ci and c 2 . For the (p,a), (p,n), (d,a), (d,n), 

(t,n), and (n,a) processes, the highest point is determined by the energy 
production at a temperature such that the corresponding most frequent 
velocity is sufficient for penetrating the Coulomb barrier. In the case 
of many reactions, the reaction cross section is larger than the scat- 
tering cross section. 

An example of the considerably less probable (p,P) reactions is the 
proton- proton reaction, indicated in the figure, with its extremely small 
cross section. 

At still higher temperatures, the duration of stay of the colliding 
particles in the hit nucleus is less and, therefore, the cross section 
again decreases so that the production curves run through the maximum 
under discussion. 

The slope of the energy production curves, on the other hand, is 
determined by the relation 

k 2 = 42.7(Z 1 Z 2 ) 2 / 3 [a i A 2 /(A 1 + A 2 ) 1 / 3 J 

that is, by the atomic number and atomic weight of the reacting particles. 
For this reason, a steeper slope is observed in figure 3 for large k 2 , 

that is, when one or both of the reaction partners has a large atomic 
number . 

This circumstance leads, for example, to the result that the energy 
rate of the hydrogen reaction chain in the colder stars at 2 to 6*10^ °K 
is proportional to T^, whereas the energy rate of the carbon- nitrogen 
reaction cycle in the hotter stars is proportional to T-*-®, as can also 
be observed in figure 3. 

At lower temperatures k 2 thus chiefly determined the energy pro- 
duction so that materials with small k 2 (like hydrogen) burn faster at 
low temperatures. 

Within a wide temperature range, the deuterium- tritium, deuterium- 
deuterium, or tritium- tritium reactions, or either of these with normal 
hydrogen, are the fastest. Heavier elements become somewhat faster only 
at extremely high temperatures. 

The numerical values shown in figure 3 hold for stoichiometric mix- 
tures of bare nuclei for which, therefore, the relation cy + c 2 =* 1 
always holds. This means that they hold for the first instant of the 
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reaction initiated in the fresh gas at the temperature T, while with 
progressive combustion naturally the cj_ and C2 decrease and, there- 
fore, the energy production per unit time decreases. 

The condition Cj_ + 03 * 1, however, is also in the first instant 

satisfied only for complete ionization of all reaction partners, but 
this is never the case at the lower plasma temperatures and thus the 
low-temperature region of the energy curves of figure 3 may be appreci- 
ably less favorable than shown. 

At the temperatures that are here of special interest above 10 7 °K 
for the light elements, the fractions of nonfully ionized particles which 
are, therefore, not capable of reaction, are generally already negligibly 
small (i.e., smaller than is desirable from considerations of the mean 
free paths of the rest mass particles) . 


Finally, it should again be recalled that figure 3 is drawn for 
fixed combustion gas densities of p « 1 gram per cubic centimeter in 
order to simplify the computations. To this density there would, for 
example, correspond for the deuteron-deuteron fresh gas of 10® °K, a 
pressure of 4.22x10® atmospheres, that is, it is necessary to reduce the 
value of the energy rates of figure 3 by about 9 powers of 10 in order 
to have technically controllable steady-state pressure relations. 


For clarifying these relations at commercially feasible, steady 
nuclear combustion conditions, the energy rate per unit volume is plotted 
against the temperature in figure 4 (according to equation (12)), for al 1 
the reactions previously considered. 

The numerical values for the unit of pressure in atmospheres are 
obtained from those of figure 3 by multiplication with the factor 



The maxima of the energy production curves per unit volume are con- 
siderably sharper than those per unit mass, also the variation between 
the different gases becomes larger, since the molecular weights enter 
quadratic ally. In spite of this, the advantage of the heavy hydrogens 
remains very clear; the reactions of the deuterium and tritium are again 
by far the fastest. The characteristic differences and difficulties of 
the technical as compared with the cosmical steady nuclear combustions 
clearly appear, which arise from the smaller pressures by about 15 powers 
of 10 and from the shorter times available for the technical combustions. 
The latter requirement makes necessary the application of higher temper- 
atures than are used by nature herself at the centers of normal stars. 
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Only at the optimal temperatures and at the same time at the highest 
possible practically attainable pressures do the half-value times of the 
hurnup and, therefore, the combustion chamber loading become of compar- 
able magnitude with those of chemical combustions. 


IV. RADIATION 

We have already established that every gas adiabatically shut off 
from its surroundings experiences thermal nuclear reactions which in- 
crease exponentially with time so that, in this case, the ignition tem- 
perature would lie at an arbitrarily low value. 

Under actual man-made conditions, the nuclear ignition temperatures 
will depend on the ratio of the energy production in the gas volume to 
the heat simultaneously given off to the surroundings. 

When this ratio permanently becomes greater than 1, the temperature 
continues to rise by itself and "ignition" occurs. 

The energy radiated can take place fundamentally both through the 
diffusion of rest mass particles and through the diffusion of photons, 
where the ratio of the geometric dimensions of the reaction plasma to 
the mean free paths of the previously mentioned types of entities is 
important. 

In this way, we come to the second fundamental difficulty of the 
practical steady nuclear combustions as compared with the -cosmic steady 
combustions. The mean free path of the particles with rest mass is 


A * 


0.177 

na Q 


(14) 


(where n a* p/kT is the number of particles per cm 3 , and Og the 

scattering cross section in cm®) . For the most favorable reaction tem- 
perature T ~ 10® °K and for the best nuclear fuel at an atmosphere com- 
bustion pressure the particle density is n » 7.13xl0-*- 3 centimeters -3 , 
which gives a mean free path A^_ * 2.5x10-*- centimeters, if we have non- 
fully ionized particles for which crg-^ ~ 10”-*-® square centimeters and 
A2 * 2.5x10® centimeters, if we have fully ionized particles where 

a S2 ~ 1°" 24 


square centimeters. 
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The corresponding number of collisions per second of a particle is 
obtained from 

8kTl Ci /A ± 

*MA 2 

which for a fresh deuterium- deuterium gas has the respective values 
Zj_ = 4.11x10® second" atmosphere"^ and Zg » 4.11xio -2 second"^- 
atmosphere"^ for the two states of ionization under consideration. 




Actually, the reaction plasma always consists of mixtures of par- 
ticles in the two states of ionization so that the actual free path A 
will lie between the two extremes Ap and Ag and are given by 

A ~ 0.177 _ 0.177 1 + n-^ng 

n l a Sl + n 2 a S2 ncr Sl 10" 8 + ri ± fn 2 

This means that for ratios n^/ng>10" 8 of nonfully ionized particles, 
the mean free path A will temporarily be proportional to ^/n-p which 

is the case for amounts of heavier and therefore nonfully ionized ele- 
ments, such as are unavoidable as natural impurities of the fuels, and 
such as occur in the purest hydrogens in consequence of the statistical 
ionization process. 


The mean free path in the reaction plasma should not, therefore, be 
expected to be larger than 


0.177xl0 16 1 + ni/n2 

n ^/ng 


(14a) 


If the geometric dimension r of the volume of the reaction plasma is 
large compared with these mean free paths , then the number of the dif- 
fusing rest mass particles within the burnup time of the plasma will be 
small. 


The product pr for the combustion gas must always be at least of 
, the order of magnitude 

- 1 + n-,/n ? 

pr •« lOOpA -2.5x10° 7-* atm cm 

n l/ n 2 


therefore, 


n-^/ng = 10" 8 , pr^ 2.5x10® atm cm 


(17) 
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The relations vill be still less favorable for the diffusing par- 
ticles without rest mass (i.e., photons). The photon cross sections of 

the nuclei are of the order of 10 - ^^ to 10”^ square centimeters, while 
the scattering cross section of photons on electrons is of the order of 

magnitude 10 “ ^ square centimeters. Thus, in a fully ionized reaction 
plasma the mean free path of the photon gas component would be many 
orders of magnitude larger than that of the gas component constituted by 
the rest mass particles. This statement is true also for the nonfully 
ionized plasma. For a gas volume whose geometric dimensions are large 
compared with the mean free path of the rest mass particles, photons 
will always diffuse and energy will be lost mainly through optical radi- 
ation, if the density of the photon component of the gas is of appreci- 
able magnitude. 


At the high plasma temperatures which are under consideration here, 
the number of the processes which lead to photon production is very high. 
Without going into their detailed description, we shall mention here 
only nuclear gamma radiation, bremsstrahlung, recombination radiation, 
radiation from transitions in nonfully ionized atoms, plasma vibrations 
and so forth, so that the loss by radiation will be practically 
continuous . 

If the dimensions of the plasma volume (e.g., of spherical shape), 
were also large compared with the free path of the photon gas , then an 
equilibrium condition would exist and the laws of black body radiation 
would hold. The energy radiated away would, in the usual way, be pro- 
portional to the area of the plasma volume and the fourth power of the 
plasma temperature. 


The energy production, on the contrary is, according to equation 
(12), proportional to the content of the plasma volume. Setting the two 
equal to each other a relation is obtained for the ignition temperature, 
or for the steady temperature, as a function of the radius r of the 
plasma sphere, if the plasma pressure p is held constant: 


6 exp(k 2 T -1 / 3 ) (^iAi) 

(koT -1 / 2 ) 2 m 2 Eci Cop 2 } 


where a » 5.68x10 ^ erg centimeter"^ degree"^ second - ^ is the Stefan- 
Boltzmann constant. 


The plasma radius has a minimum at the point where, in figure 4, 
the tangent of the T^ law touches the energy rate curves, which for 
deuteron-deuteron reaction is at about 10 7 °K. The magnitude of this 
minimum radius in this case and at a plasma pressure of 100 atmospheres 
17 

is about r +* 10 centimeters, therefore, of astronomical order of 
magnitude . 


ATT* 
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For a particle density of n ~ 10^ centimeters"'-’ and an optical 

cross section of a ~ 10" zz square centimeters, the mean free path of 
the photon gas becomes A = 0.177/n a or the thickness of an optically 
infinite layer of the plasma is r® * 4.6/ncr which is of the order of 

10 17 centimeters. Therefore, in stationary plasmas of practical dimen- 
sions, no black radiation can occur in spite of the continuous radiation 
because the optically infinite layer thicknesses required for this pur- 
pose cannot be attained. 

It will nevertheless, in this case, be necessary to consider approx- 
imately gray radiation whose intensity ratio to black radiation is given 
by the relation 


l/l 0 « 1 - e -rn0 ~ rna 


(19) 


The dependence of the effective optical cross section a on p 
and T in the absorption constant no we shall assume, in the absence 
of better data, to be the expression given by Kramers, Eddington, and 
Gaust (ref. 2) for the plasma at the interior of the sun: 


a ~ 6. 60(l + c H )(l - eg - cue) “ 


*i z iAi 


Zci/Ai 


p0.7b rp-3,5 


6.60(1 + c H )(l - c H - c He ) - — — 


J c iZ|/Ai 


( 20 ) 


(k/Mx 2 c i /A i ) 


1.75 


pi .75 ip— 1.25 


so that the radiation intensity of the plasma is 


I ~ iQrna = 


3.97xl0 Z4 ar(l + Cg)(l - eg - cg e ) 


^± 4 /h 

i 

(k/MxZc./Ai ) 1 * 75 

i 


a. 75 




( 21 ) 


At constant plasma pressure we thus arrive at the surprising result 
that the radiation intensity temporarily drops with increasing tempera- 
ture even somewhat more sharply than the plasma density, until at a 
higher temperature, with the onset of intense bremsstrahlung, mass radi- 
ation, and finally the vanishing of the boundaries between matter and 
energy, the radiation intensity again rises. 
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From equations (12) and (21) is now obtained the very important 
ratio of the energy radiation loss to the energy production for a plasma 
sphere of radius r: 


i = 


4r 2 rtl 


4/3x r 3 rtAE" 


= 3nal 0 /AE" 


1.192X10 25 a(l + c H )(l - c H - c He ) 


'H 


Ec l c 2 k l k 2 


(l^i/Ai) 0 * 25 ^ c iZi/AiXp' 0 - 25 T 1 - 417 expC^T -1 / 3 ) 


( 22 ) 


which is independent of the dimension- r of the reaction plasma, depends 
little on the plasma pressure p and is essentially determined by the 
temperature T. 


The ratio t, for constant pressure p, therefore, has a minimum 
with respect to T at 

T = (~^—f (22a 

\4.251 J 

a value which is again independent of the pressure and remains somewhat 
above the optimum temperature of the energy production as shown by 
equation (l2a) . 


The magnitude of the minimum ratio £ is 


3 min 


1.77X10 24 a(l + c H )(l- c H - c He ) ^ ^ 2/> /'kV 

, .-2.251 C i z i/ A i( mZ/ 

1 c 2 K 1 K 2 1 V i 



p-0’25 


(22b) 


The Kramers relation (eq. 20) assumes, perhaps somewhat optimis- 
tically, that plasmas which consist of only fully ionized particles like 
hydrogen and helium do not radiate, so that the actual radiation is 
caused by the unavoidable impurities of the nuclear fuels with heavier 

elements, or to the 10 - ^ fraction of nonfully ionized particles as shown 
by equation (16) , to limit the mean free path length. 


In this way, the factor (l - Cg - cg e ) in the preceding equations, 

becomes equal to 10 -3 . With this somewhat arbitrarily chosen number, 
the radiation ratio £ is plotted in figure 5 as given by equation (22); 
also given are the optical cross section by equation (20) , and the abso- 
lute radiation intensity I by equation (2l) at the deuterium- deuterium 
plasma pressure of 100 atmospheres, where for the computation of I, the 
radius of the plasma sphere was assumed to be r » 1 meter. 
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The plot of the number of particles n per cubic centimeter and 
per atmospheric pressure against the temperature shows that for nuclear 
reaction plasmas, the densities considered are such as would otherwise 
he denoted as a high vacuum. 

The plot of the optical cross sections a (referred to the relative 
mass concentration of the particles which are neither hydrogen nor helium) 
shows the extremely strong drop of this average optical cross section 
with increasing temperature and its relatively small dependence on the 
pressure. In the region of the most favorable reaction temperatures of 
good nuclear fuels a lies at about 10 -33 square centimeters. Corre- 
spondingly, the optical absorption constant nd likewise drops approxi- 
mately with t - 5*25 an £ resu i-ts in the decrease of the radiation inten- 
sities I decrease approximately as t - -'-*^ when the optical layer 
densities rncr become less than 4.6. For rad >4.6, the intensity is 
practically the same as that for black radiation, also plotted as the 
asymptote to all the radiation curves. 

Finally, the most important plot of the radiation ratio £ shows 
that under steady, realizable conditions at temperatures below about 

10 7 °K the radiation always remains larger than the energy production, 
so that below this temperature a steady nuclear combustion cannot be 
obtained. At higher temperatures, the energy production may become 
larger than the radiation if the proportion of heavy or nonfully ionized 
particles in the plasma is sufficiently small, of the order of magnitude 

lCT 3 to 10" 7 . 

The other curves in figure 5, for various values of p, r, and 
(1 - c H - c He ) , require no detailed discussion but serve to complete 
essentially this section on the radiation relations of steady technical 
nuclear combustion plasmas. 

These fundamental considerations show that with decreasing content 
in the plasma of the nonfully ionized particles the radiation of photons 
rapidly decreases; also, the radiation from the massive particles in- 
creases so that a minimum of the energy discharged, with respect to the 
content in nonfully ionized particles, is obtained. 

In this way, the conditions under which it may be possible to avoid 
having the energy production of the nuclear combustion be initially lost 
to the surroundings, instead of being used for the preparation of the 
fresh gas, are approximately outlined, which evidently, is an essential 
condition for steady nuclear combustion. 
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V. MIXTURE PREPARATION 

The mixture preparation in steady nuclear combustions embraces the 
same phases as those in chemical combustions (ref. 3) , in particular, 
the fuel injection, the jet atomization, the droplet formation, the drop- 
let evaporation, the macromixing, the micromixing, the heating up to igni- 
tion, and, finally, the ignition itself. 

We confine ourselves here, owing to its extraordinary energy takeup, 
to the particularly important partial process of the heating up to the 
ignition, which includes also the total ionization processes of the 
fresh gases, since in the case of nuclear reactions, the ignition tem- 
perature lies above the ionization temperature. 

As previously mentioned, chemical and nuclear chemical reactions 
take place continuously in every gas of our environment and consequently 
in the combustion gases of our jet engines although this occurrence is 
so rare it is generally unobserved and is of little importance from the 
energy standpoint. The great majority of the collisions of the gas par- 
ticles, due to their thermal motion, concerns only the electron envelopes 
of the colliding atoms without any occurrence of direct collisions be- 
tween nuclei. 

This condition of thermal nuclear reactions is satisfied, to a 
great extent, only when combustion gas temperatures are so high that the 
greater number of the gas atoms is fully ionized and are completely 
freed of their electron envelopes. Investigations at these high com- 
bustion temperatures, particularly of electrons, must, strictly speaking, 
include the methods of relativistic thermodynamics, but we shall be sat- 
isfied here with several approximate considerations on the basis of 
classical thermodynamics. 

To have an equilibrium of the ionization processes requires an equi- 
librium between the collision ionizations occurring in the first instant 
due to collision between nonionized atoms and the ions and electrons 
taking part in further collision ionizations. There must also be a 
regeneration of the high velocity end of the Maxwell distribution that 
was consumed in the first stages of ionization. Accompanying the in- 
creased ionization there is increased recombination, thereby acting to 
reduce the degree of ionization. Equilibria between all of these 
processes require about 10 9 collisions. Thus, equilibrium will not be 
set up during the burnup of fast nuclear chemical combustions as experi- 
enced in the case for chemical combustions. 

During fast nuclear combustions, the pressure- dependent ionization 
equilibrium (such as is attempted to be described by the Eggert-Saha 
equation) does not have to be considered, but is actually maintained in 
the slowly burning plasmas of the cosmic nuclear combustions. During a 
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considerable part of the burnup, we shall mainly have to consider only 
the pressure- independent primary spontaneous ionization which, depending 
on the combustion gas pressure under consideration, can lead to smaller, 
equal, or the same degrees of ionization as the pressure- dependent 
ionization equilibrium. We shall next consider the primary ionization 
in more detail. 

Collision ionization is to be expected in the case of thermal col- 
lisions between particles of different mass A-^M and A2M if the energy 

AfAg/CAi + Ag)><Mw^/2 = AMw 2 /2, which vanishes in the inelastic collision 

in the collision system, is equal to the ionization energy I. 

The average required collision velodity w^ required for ionization 
is 


w ± « V2I/AM (23) 

For equal masses of the colliding particles, this means that the 
relative kinetic energy of a particle A^Mw^/2 must, on the average, be 

at least twice as large as the ionization energy which vanishes in the 
collision system so that the remaining half of the kinetic initial energy 
is still able to satisfy the requirements of conservation of momentum. 

Within the Maxwell-Boltzmann distribution of the thermal particle 
velocities it is thus necessary to have the most frequent velocity 

c = ^/2kT( 2c. /A. ) /M (24) 

1 X 


attain at least the required collision velocity w^, and from this the 
required temperature of the gas for full ionization is obtained as 


kAZCi/Ai 


(25) 


The energy I required for the complete ionization of an atom de- 
pends on the atomic number Z, approximately, according to the relation 


I = 2.18X10- 11 z 2 - 418 [erg] 


(26) 


Substituting this in equation (25) there is finally obtained 


2.18X10" 11 Z c,Z?' 418 


11 


T a 


kA Zc i /A i 


(27) 
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Below a gas temperature of 5x10^ °K and from consideration of ion- 
ization there is, therefore, no expectation that any gas or any type of 
hydrogen would produce a reaction yield of any technical importance when 
dealing with a fast, short time combustion process. 

At the temperature of the centers of normal stars ( from 2° to 6° * 10 ' ) 
only the first eight elements of the periodic system up to oxygen 
can be expected to give thermal nuclear reactions, a choice which holds, 
in quite the same way for the chemical reactions in chemical rockets. 

Below the full ionization temperatures represented by equation (27), 
only small fractions of the gas mass can take part in the nuclear reac- 
tions. According to the Maxwell distribution curve, the number of reac- 
tions become smaller when further removed from the actual gas temperatures 
below the values of equation (27). 

The numerical fraction of the gas particles N*/N, which are faster 
than Wj_, from which, therefore, complete spontaneous ionization is to 
be expected, is approximately 



In accordance with the approximate nonrelativistic treatment, 

N'/N «s 1 only for c * «°, that is, for infinitely high temperature. 

For Wj^/c * 1 there is found N'/N » 0.58, that is, at the temper- 
atures represented by equation (27), a majority of 58 percent of the 
particles always attain the collision velocity required for complete 
ionization, at half the temperature the number is still 25 percent, at 
T/4 the number is 4 percent, and at T/8 considerably below 1 percent. 

The fraction of completely ionized particles of N'/N * (l - lO"^) 
in figure 5 is, according to equation (2o,, attained for a temperature 
which is about 25 times as great as that according to equation (27). 

Only at these extremely high temperatures, especially for the heavier 
elements, does the amount of complete ionization rise and the radiation 
therefore drop to an extent that makes continuous nuclear combustions 
possible. 

The next problem concerns the energies required to heat the fresh 
gas to these high combustion temperatures. 
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We first assume that the heating process is adiabatic, t h at is, 
neglect energy losses to the outside due to diffusing photons or elec- 
trons or due to convective heat transfers to the boundary walls by atomic, 
ionic, nuclear, electron, or photon collisions against the wall. 

This assumption we are justified in making under certain conditions 
according to the preceding section. 

This assumption will again closely apply for very rapid heating 
processes for which the partial pressures of the previously mentioned 
five most essential combustion gas components also do not have time for 
establishing their own equilibria. 

In spite of these simplifying assumptions, the problem of the deter- 
mination of the specific heats of the combustion gas is also of a static 
character although, because of the predominating effects of the level 
excitations and their decay, the specific heats become functions of the 
time. 


The Boltzmann equipartition law, as a consequence of the classical 
Boltzmann statistics, still possesses within the framework of the quantum 
theory, the significance of a limiting law which is valid for sufficiently 
high temperatures , as under the conditions we are considering. It does 
not hold, however, for the photon gas which leads to the Rayleigh- Jeans 
radiation law, valid only for long wavelengths. 

The Boltzmann equipartition law states that each degree of freedom 
of a mechanical system, whose energy depends quadratieally on the veloc- 
ity or coordinate has, on the average, an associated energy kT/2. The 
mean energy of a gas atom, with its three spatial translational degrees 
of freedom, therefore, amounts to 3kT/2. 

If an atom of this energy is supplied with the ionization energy 
Ip, an electron is split off which, in accordance with experimental re- 
sults, after only a few collisions, comes into translational energy 
equilibrium with the remaining gas, that is, adjusts to the Maxwell- 
Boltzmann velocity distribution and, therefore, in turn takes up the 
mean energy content of 3kT/2 so that, referred to the original atom, 
a total energy of (3kT/2 + Ip + 3kT/2) must have been made available. 

If the atom ionizes multiply, for example Z times, the total mean 
energy increases to 

Z Z 

U z * 3 kT/2 + £ I + 3ZkT/2 « £ I + (Z + l) 3kT/2 
1 .1 


( 29 ) 
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The specific heat Cp at constant pressure and per unit mass 
approaches asymptotically, therefore, with infinitely increasing temper- 
ature, the limiting value 


C P 



z i + 1 5k 
Aj_ 2M 


(30) 


while for very low temperatures it corresponds to the value 


c 


P 



5k 
Ai 2M 


(31) 


The transition curve between these two values depends on the rapid- 
ity of the heating. For the short time heating here assumed, the trans- 
ition curve can be given with the aid of equation (28) for intermediate 
temperatures. For medium temperatures, there is generally obtained a 
weak maximum, due to the ionization energy of the mean specific heats. 
The maximum may easily lie far above the limiting value as shown by 
equation (30) . 


While the initial values of the specific heats per unit mass thus 
decrease mo no tonic ally and hyperbolic ally to zero with increasing atomic 
weight, the final values at high temperatures approach somewhat less 
regularly the final value Cp ~ 2 calorie gram” 1 degree" 1 , since the 

value (Z + l)/A drops toward 0.4 with increasing atomic weight. 

The largest changes in the specific heats occur for the lightest 
elements, especially the various hydrogens, as shown in figure 6. 

If we are not considering short time processes, but assume that a 
thermal equilibrium is established, then the horizontal branches of the 
curves of the mean specific heats remain unchanged. Only the S- shaped 
transition portion, which represents the ionization processes, shifts 
in a manner that depends on the temperature. At very high temperatures, 
we also have, in this case, the same enthalpy as in short-time processes 
if the level excitations and reradiations are not considered. 


With these assumptions and with the aid of figure 6, we can con- 
struct an enthalpy- temperature diagram, as shown in figure 7 for example, 
for stoichiometric fresh gas mixtures of hydrogen and deuterium with the 
first elements of the periodic system. 

The enthalpies, up to the corresponding value of the heat of the 
nuclear reaction, are drawn as a solid line; beyond that the line is 
dotted. In addition, the optimum reaction temperatures shown in column 
12 of table I are plotted. 
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The heats of reaction for the hydrogens are quite sufficient for 
heating to the most favorable reaction temperatures, but for the reac- 
tions of the heavier elements, this is not always the case. Since it 
is also for the heavy elements that the radiation during the reaction is 
considerably higher than for the hydrogens, the fundamental possibility 
of continuous nuclear combustion for these heavier elements is very 
questionable. 


VI. APPLICATION TO JET PROPULSION 


From the preceding investigations, which undoubtedly are associated 
with a number of uncertainties of considerable importance, the following 
most important conclusions for the achievement of continuous nuclear 
combustions can be drawn: 

(1) The fuel of primary consideration appears to be deuterium of 

the highest purity, with respect to heavier elements. The reasons for 
these choices are: the small half-value time of burnup, high energy 

production per second per unit volume of the combustion chamber, small 
radiation intensity of the reaction plasma and high heat of reaction in 
comparison with the enthalpy of the fresh gas at the optimum reaction 
temperature and economical operation. 

(2) The combustion chamber loading increases as the square of the 
pressure employed, but even for combustion pressures of 100 atmospheres^ 
attains only values of the order of magnitude of 10 4 kilocalories meter"3 
second - -*-, and, therefore, is less than that of chemical high pressure 
combustions . 


(3) The minimum dimensions, even of high-pressure combustions, are 
of the order of many meters since the mean free paths of the high- 
pressure reaction plasma are of the order of centimeters. 


(4) For a cooling of the burned up plasma considerable below the 
most favorable combustion temperature of 10® °K there occurs considerable 
recombination and radiation of extraordinary intensity which, at about 
10 5 °K, may cause the optical radiation to amount to a considerable part 
of the total energy production. 


From these circumstances it is evident that nuclear combustions for 
stationary units or those used in ships are more readily realizable than 
those employed in aircraft or spacecraft, where in the former, the min- 
imum outputs required are about 200,000 horsepower (several milligrams 
of fuel consumption per second) , which can more easily be exceeded. 
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For jet propulsion units, three fundamental processes for utilizing 
the combustion gas enthalpy of lO 11 -* kilocalories per kilogram have been 
considered: 

(l) The direct release of the combustion gas such as that used in 
the familiar chemical rockets (i.e., the so-called pure atomic rockets). 
The exhaust velocity would be of the order of magnitude 10^ meters per 
second, and the specific fuel consumption about 10“^ kilogram per 
second,' so that a combustible layer of protective gas or protective 
liquids for the furnace and nozzle walls appears tolerable. The ignition 
of the reaction plasma with the cooling layer, which would occur if the 
evaporation proceeded more rapidly than recombination, could be restricted 
to a possible mixing process of the plasma with the protective gas. 

This circumstance could make the pure' atomic rockets appear in a new 
and not altogether hopeless light. The radiation of the reaction plasma 
itself which, as shown in figure 5 is 1000 calories per square centi- 
meter second, could be further controlled in the combustion chamber 
walls except when the penetrating x-ray radiation predominates. 


The convection product Hpc that determines the convective heat 
transfer and has a value of about 4x10^ kilcalories per square meter 
is, in spite of the small gas densities, about 10^ times larger than that 
used in chemical high-pressure rockets and does not permit any antic- 
pation of controlling the heat transfers at an unprotected wall by simple 
wall cooling. On the other hand, because of the unusual kinematic vis- 
cosity ca/ 3 of the plasma flow and its correspondingly low Reynolds 
number r/A ~ 10^, there are favorable prospects for the mixing zones of 
protective gas and plasma which suggest a compromise solution between 
pure and thermal atomic rockets as considered by H. J. Kaeppeler (ref. 4). 


(2) The next method for utilizing the nuclear plasma enthalpy con- 
sists in completely mixing the hot plasma with the surrounding air 
(atomic air jet units) or with other inert working gases (thermal atomic 
rockets) , thereby heating these gases to a temperature lying only one 
order of magnitude above that of chemical propulsion units. Thermal 
atomic rockets of this kind evidently differ only in degree from the 
previously mentioned pure atomic rockets having a protective gas layer, 
although for both admixture applications, it is preferable to have the 
lowest plasma pressures at which the energy production is just sufficient. 
Moreover, in both cases during the mixing, the temperature range of 
extremely high radiation intensity must be run through so that no gen- 
eral conclusions can be drawn without a detailed study of these special 
relations, which have been taken into consideration by I. Sanger-Bredt 
and H. J. Kaeppeler (ref. 4) . 


In any case, it may be said that the energy radiation cannot become 
larger than the energy production and the value for a combustion pressure 
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of 10 atmospheres and for a radius r(m) of a spherical chamber would be 

about I = 23 r kilocalories meter -2 second -1 , so that, even for a 10 
meter combustion chamber radius, the radiation is less than 0.023 kilo- 


calories centimeter - ^ second - - 1 , which is feasible to handle and may per- 
haps be of interest for stationary nuclear combustion units. 


(3) Finally, a propulsive method utilizing the high radiation 
intensities that were the cause of the cooling reaction plasma in the 
other two methods will be considered (the so-called photon rocket) 
(ref. 5). 


If this cooling of the reaction plasma occurs in a closed space 
with, for example, central reaction zone and surrounding cool gas jacket, 
that is, without appreciable gas flow, it should be easily possible to 
arrange for the greatest part of the plasma enthalpy to be converted 
into photon gas and radiated. 


The radiation intensity again depends linearly on the lamp radius 
r(m) and on the square of the combustion pressure p(atm) and for the 
triton- deuteron reactions is supposedly 

I * 3.7rp^, kcal m -2 sec -1 (32) 

corresponding to a radiation pressure of 

p s * 5.6xl0 -1 ^ rp^, atm (33) 

In order, therefore, to attain radiation pressures of an order of 
magnitude that is of technical interest (i.e., at least 0.01 atm), the 
combustion pressures for a lamp radius of r «■ 1 meter must amount to 
p * 4.2*10^ atmospheres. 


For photon rockets, therefore, in the same way as for atomic air 
jet propulsion units, it will still be necessary to seek energy sources 
of higher reaction velocity in order to reduce the technical difficulties. 


REFERENCES 

1. Gamow, G. , and Critchfield, C. L.: Atomic Nucleus and Nuclear Energy 

Sources. Clarendon Press (Oxford), 1950. 

2. Gamow, G., and Critchfield, C. L. : 1. c., p. 270. 

3. Sanger, E. : Brennstoff-Chem. , 32, 1, 33, 1951; 35, 41, 1954. 



32 


NACA TM 1405 


4. Sanger- Bredt, I.: Probleme der Weltraumforschung. Biel, Laubscher, 

1955 . 


Kaeppeler, H. J. : Probleme der Weltraumforschung. Biel, Laubscher, 

1955. 


5. Singer, E. : Probleme der Weltraumforschung. Beil, Laubscher, 1955. 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics 


TABLE r. - CHARACTERISTIC VALUES FOR STOICHIOMETRIC TWO-PARTICLE NUCLEAR COMBUSTION OF CHARGED LIGHT PARTICLES 


NACA TM 1405 


33 


Maximum energy rate per unit 
volume at constant plasma 
pressure 

< 
i— ( 

CO 1 

Pi O 

0) 

X w 
■? 

3 § 

i — i 

a 

o 

CO 

CO 

1 

o 

rH 

X 

8 

m 

i 

O 

rH 

X 

LP 

CO 

CO 

f- 

3 

5 

r~ 

to 

00 

1 

o 

i—t 

5 

LQ 

rH 

ip 

i 

O 

H 

ft 

rH 

IP 

1 

o 

i — i 

8 

CD 

rH 

LO 

( 

3 

ft 

Tii 

CO 

O 

rH 

o 

rH 

ft 

ip 

ip 

LP 

i 

O 

X 

00 

r- 

ip 

T# 

o 

rH 

X 

to 

00 

rH 

•HI 

1 

o 

rH 

8 

to 

to 

co 

1 

o 

rH 

X 

8 

rH 

Til 

1 

o 

H 

X 

t" 

CD 

CO 

tIi 

o 

H 

X 

8 

CO 

to 

0 

rH 

1 
tH 

3 

Pi Ut 

V°° 

r~ 

O 

rH 

X 

rH 

to 

r- 

% 

rH 

X 

CO 

o 

rH 

% 

H 

ft 

H 

rH 

b 

rH 

6 

to 

ft 

CD 

o 

rH 

ft 

H 

rH 

CD 

o 

rH 

ft 

CO 

f-H 

s 

$ 

rH 

CO 

b 

rH 

to 

b 

rH 

X 

CO 

ip 

to 

CD 

o 

rH 

P- 

% 

E 

LO 

rH 

b 

H 

X 

to 

00 

rH 

CD 

o 

8 

o 

CO 

CD 

3 

X 

3 

CO 

b 

rH 

3 

to 

CO 

3 

Maximum energy rate per unit mass 
at constant plasma pressure 

s 

a> *n -» -p 

*iH 3 a) +jco 
H H 1 Pi 

as a -h ;Oo 

S > +j EH ti 

IQ 

to 

CO 

o 

i — i 

X 

CO 

to 

o 

rH 

ft 

CO 

CD 

b 

i 

CD 

IP 

'o 

rH 

X 

to 

CO 

CO 

rH 

1 

o 

rH 

ft 

ip 

tii 

% 

8 

rH 

IP 

rH 

3 

V 

K 

CD 

rH 

b 

rH 

X 

CO 

CO 

rH 

rH 

o 

X 

P 

rH 

CO 

o 

ft 

rH 

rH 

rH 

o i 

3 

to 

to 

rH 

1 

o ! 

rH 

ft 

o 

H 

rH 

1 

CO 

r- 

ft 

H 

o 

E 

LO 

rH 

rH 

o 

rH 

8 

CO 

ft 

rH 

rH 

1 — \ 
a 

aj 

Pi W 

X H 

o 

S 

-p 

CtJ 

T* 

rH 

1 

o 

rH 

X 

in 

CD 

CO 

3 

s 

CO 

r- 

O 

rH 

ft 

to 

rH 

b 

rH 

8 

ip 

to 

o 

b 

rH 

ft 

cn 

rH 

CD 

o 

rH 

ft 

to 

to 

% 

E 

r- 

rH 

b 

ft 

CO 

to 

1 

1 

to 

rH 

H 

O 

X 

r~ 

rH 

Tf 

o 

b 

rH 

ft 

LP 

to 

rH 

ft 

co 

ip 

H 

b 

H 

ft 

LP 

CO 

H 

ft 

CD 

rH 

rH 

b 

X 

CO 

3 

& « 
- o o 

EH 

b 

rH 

X 

3 

to 

% 

rH 

X 

to i 

CO 

t# 

b 

rH 

ft 

ft 

o 

rH 

X 

CD 

CO 

rH 

b 

rH 

8 

co 

ft 

CD 

o 

H 

IP 

°b 

H 

8 

rH 

cn 

o 

b 

s? 

ip 

Til 

rH 

O 

b 

X 

rH 

o 

b 

H 

8 

CD 

CO 

% 

E 

to 

to 

b 

rH 

ft 

ip 

r- 

b 1 

X 

to 

ip 

co 

b 

8 

n- 

co 

% 

rH 

X 

LP 

3 

CD 

CD 

Maximum energy rate per 
unit mass at constant 
plasma density 

rH 

•* V 

CL (1) 

« 
X , 

d i-H 

aJ 

a 

tO 

l 

e 

M 

r- 

O 

rH 

X 

r- 

CO 

b 

rH 

ft 

t- 

ip 

rH 

3 

5 

CO 

H 

O 

ft 

t* 

rH 

to 

rH 

o 

rH 

X 

IP 

Til 

co 

f- 

b 

H 

ft 

to 

CO 

r- 

b 

(H 

8 

Til 

00 

3 

O 

H 

ft 

CD 

to 

8 

to 

to 

CD 

b 

H 

§ 

to 

r- 

b 

rH 

ft 

CD 

CO 

co 

b 

rH 

8 

rH 

cn 

Oi 

b 

rH 

g 

ft 

E 

CD 

00 

00 

b 

rH 

ft 

co 

ft 

00 

Pi W 
O O 
EH 

b 

rH 

X 

CO 

CD 

ft 

rH 

§ 

ip 

CO 

% 

rH 

ft 

to 

f- 

1 

X 

to 

CD 

CO 

o 

H 

o 

rH 

X 

CO 

ip 

r» 

s 

o 

H 

ft 

r- 

r- 

rH 

| 

rH 

rH 

•4 

8 

CO 

CO 

-o 

8 

CO 

co 

=i 

8 

IP 

T^ 

•8> 

8 

00 

CD 

0 

1 
s 

H 

H 

ft 

lO 

rH 

E 

ro 

rH 

o 

b 

rH 

g 

H 

r- 

Material constant of 
Gamow equation 


to 
CO i— 1 

X M 

■3 

b 

rH 

X 

CD 

to 

to 

■a 

E 

r- 

to 

b 

1 — l 

ft 

CO 

to 

b 

rH 

ft 

rH 

IQ 

b 

s 

CO 

% 

8 

ip 

co 

1 

rH 

t 

CVJ 

rH 

1 

CO 

CO 

rH 

t> 

rH 

ft 

LO 

rH 

8 

CO 

T^ 

■s 

g 

ft 

b 

H 

g 

rH 

3 

g 

rH 

! 

CD 

Tii 

to 


r — 1 

o 

0) 

to 

hO 

u 

b 

rH 

2 

CO 

°b 

rH 

X 

co 

Tt< 

CO 

r- 

°b 

H 

6 

co 

CO 

* 

r- 

b 

H 

ft 

CO 

rH 

rH 

b 

rH 

8 

r- 

00 

o 

*°o 

rH 

8 

to 

CO 

co 

b 

rH 

ft 

o 

CO 

r- 

1. 

g 

to 

CO 

CO 

% 

t 

to 

CO 

to 

■& 

S 

Til 

rH 

b 

rH 

ft 

co 

T# 

CO 

■5 

o 

rH 

t 

co 

b 

rH 

ft 

to 

rH 

* 

CO 

b 

rH 

ft 

to 

rH 

< 

rH 

b 

H 

ft 

r- 

t 

LP 

Resonance 

width, 

r, 

ev 

CD 

H 

o 

rH 

X 

CO 

rH 

o' 

rH 

o 

rH 

b 

rH 

ft 

s 

4 

b 

rH 

3 


P“ 

o 

H 

3 

8 

b 

H 

b 

rH 

b 

H 

b 

rH 

Til 

Heat of 
reaction, 
E, 

Mev 

(cal g _1 ) 

6"" 

rH 

o 

CO rH 

rH <D 
rH 

o' 

H o 

rH 

X 

"i* 

ip co 
to 

iH 

b 

H 
X 
IP 
O H 
co * 

rH 

rH 

b 

<J) H 
H • 
H 

o' 

H 
O 
rH 
00 X 
• ip 
to CO 

rH 

o' 

t 

CD 
rH • 
tJI 

CD 

o 

rH 

^ft 

CO O 
IP 

s' 

o 

rH 

oo oo 

A 

o' 

rH 

o 

rH 

00 X 

■ CD 

00 to 

rH 

°o 

Til rH 

®ft 

Til CD 
CO 

cT 

H o 

H 
to co 

ft 

o' 

b 

rH 

3ft 

CD 

o' 

b 

rH 

T!ft 

co to 
co • 
to 

o' 

b 

H 

ft 

H 3 

to 

o' 

rH 
O 
H 
X 
CO 
H CO 
rH » 
Til 

to 

React ion 

% 

a; 

oi 

rH 

m 

to 

>~ 

pi 

CO 

w 

t* 

t* 

jS 

+" 

<D 

kT" 

w 

to 

M 

ffi 

? 

pi 

to 

•H 

1-3 

Til 

ffi 

CO 

b 

vti 

r~ 

■»H 

cO 

■H 

Ml 

"d 

cn 

£ 

rH 

rH 

O 

5- 

3 

o' 

i-H 

PQ 

Til 

s 

CO 

a 

S' 

PQ 

CO 

rH 

CO 

"d 

pi 

iT 

& 

to 

il 

b 

Tl* 

JS 

'a 

ft 

w 

Tt* 

«d 

to 

■H 

•-3 

LP 

jS 

"d 

nq 

h- 

•H 

•J 

t* 

""fl 

CO 

-p 

b 

CO 

p 

o> 

I 

rH 

CO 

to 


IP 

to 

r- 

00 

CD 

3 

rH 

rH 

3 

to 

H 

rH 

in 

rH 

rH 



MCA TM 1405 



W 


Figure 1. - Maxwell- Boltzmann distri- 
bution dltf/dw of particle ve- 
locities in nuclear combustion 
plasma; variation of cross section 
a of thermal nuclear reactions 
with particle velocity w and 
curve of reaction frequency, 
a dN/dw. 
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Figure 6. - Mean specific heats of ionizing gases 
in short time processes. 






